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ABSTMCT 


The present work has been undertaken to devise 
methodologies for determining the optimal locations of new 
facilities within a region all over which demand for these 
facilities exists# Models for the location of single new 
facility and the s3.multsm.eous locat3.on of multiple new 
facilities with interaction sunong the new facilities have 
been developed* 

The thesis, for the most part, centers arotmd the 
facility location models with continuous demand, but a proba- 
bilistic version of the discrete demand single facility 
location problem has mlso been developed because this 
probabilistic version is very much s3jnilar to the continuous 
demand single facility location problem as far as its nature 
Sind method of solution are concerned. In this model, the 
co-ordinates of the destination points are assumed to be 
random variables instead of deterministic as assumed in ttie 
previous works in location theory. 

All the above stated models have been solved for 
three types of distances, viz. rectilinear, squared Euclidean, 
and Euclidean, 
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Chapter 1 


IlTTBOroCTIOH 

Selection of most profitable locations for plants 
or other fac3.1ities has been a subject of analysis for 
centuries. It has attracted a wide diversity of practitioners; 
including applied mathematicians interested in the mathema- 
tical characterization of proposed algorithmic solution 
mod-els, geographers interested ixi the spatial patterns of 
mass human interactions 5 urban and regional economists attenw 
pting to deal with problems of city land use and agricultural 
location? and planners, systems engineers, and operations 
researchers concerned with solving specific problems relating 
to the location of a particular facility, Economists, uiban 
planners, management scientists, architects, regional scien- 
tists, home economists, and engineers from several disciplines 
have discovered a coromonality of interest in their concern 
for the optimal location of facilities, 

Sven though facility location problems have received 
considerable interest over the years, it was not until the 
emergence of the techniques of operations research that the 
subject was studied in a quantitative manner. In fact, 
before 1960, the approaches to tackle these problems relied 
heavily on intuition and engineering judgement, and used 
only qualitative methods giving no or very less consideration 
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to the possib5JLity of discovering analytical solution 
procedures* 

The past ten or fifteen years have seen a signifi*« 
cant development in what has come to be known as modern 
location theory to provide the facility analysts with new 
techniques, approaches and philosophies for the solution of 
facility location problems and to take a more quantitative 
approach than was commonly taken before. In fact, there are 
quantitative aspects of the facility location problems that 
cannot be reckoned accurately through intuition alone* 

Modem location theory provides an in-depth treat- 
ment of a relatively limited number of aspects of the locational 
problems as compared with the traditional approach, which 
examines a broad rai^e of problems in a less scientific 
fashion. It deals with a smaller set of problems because it 
is still in its early stage of development. Many theoretical 
and applied problems remain to be solved, 

Sven though heavy emphasis is given, in modern 
location theory, to analytical approaches, one should rot 
overlook tl^ importance of qualitative aspects of locational 
problems* It should be realized that the analytical approach 
yields the solution to the model but not necessarily to the 
problem* The solutions obtained through analysis serve only 
as aids in decision making. There remains a number of 
non-qiiantifiable f actors j whi<^ aust be:(Cons5-dered alongwith 


scheme is presented in Chapter 2, It is bel5-eved that the 
proposed classification scheme will provide a modest push 
to the need of a bird* s eye view of the location theory 
problems. 

Although considerable amount of work has been done 
by previous researchers on the problem of locating the 
facilities when there are no constraints upon the area in 
which they can be located, yet the models given by them 
cannot be used when there are a very large number of d^and 
points needing the services of the new facilities. An attempt 
has been made in this work to develop mathematical models 
and solution procedures for continuous demand location 
problems* 

In Chapter 3 we shall describe the problem in which 
only one new facility is to be located and in Chapter 4 we 
shall deal with the problem of simultaneous location of n 
new facilities* In Chapter 5, the probabilistic version of 
the single facility location problem will be discussed. 

Our approach will be to first give a brief account of what 
is already available in literature concerning these models 
and then to describe in detail the proposed extension of the 
models* ■ 


Chapter 2 


CLASSIFICATION OF F ACILITY L OGAT IOH PIPBLBMS 

The development of literature on facility location 
problems is somewhat haphazard. The available literature on 
location theory deals with a large number of location problems 
and there seems to be little co-ordination bet^-reen tliem'y 
making it difficult to find where a particular problem stands 
in the field of location theory, A thorough review of the 
exist3.ng literature on location theory indf.cates that the 
location problems can be grouped by using the following 
scheme : 

1, Classification based on the solution space 
available, 

2, Classification based on the optimization criterion, 

3, Classification based on the number and character- 
istics of the facilities to be located, 

and 4, Classification based on the planning horizon. 

The various types of problems which belong to four 
classifications listed above are discussed in the following 
paragraphs, 

. Classification Based on the Sol u tion Space Avai lable j 

This criterion classifies the facility location 
problems (ILP*s)‘ according to the type of area which is 
available to us for locatir^ the new facilities. There are 
three categories of problems which stem from this 


classification. The three categories are : 

a) Continuous space location prohlems, h) Discrete 
space location problems j and c) Location problems on graphs. 
Continuous Space Lo cation Prob l ems - Facility 
location problems in which new facilities can be located 
anywhere in the plane come under this category* The location 
of small depots, schools, shops and all those facilities which 
do not require a large floor area or some special commodity 
which is not available everywhere are some of the typical 
examples of continuous space location problems* 

For locating the facilities in continuous space, 
one has to deteimine the distance betvreen various facilities 
for evaluating the total materials handling cost. Depending 
on the requirements of the problem, distances, are measured 
n.n several ways. Most commonly used methods are : 

(i) Euclidean, (ii) Rectilinear, and (iii) Squared 
Euclidean, 

(i) Euclidean Dist ance. The most obvious distance 
between the facilities is the straight line, or Euclidean, 
distance. If the co-ordinates of ti-zo facilities are 
(x;j_, y3_) and (x2j y2) ^spectively, the Euclidean distance 
betvreen them is defined as 

■ 1 ■ 

d= 

Euclidean distance is used for some netirork loca- 
tion problems as well as some instances Involving conveyors 
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and air travel. Some electrical "v/iring problems and pipe- 
line design protleras are also examples of Euclidean distance 
problems* 


(ii) R ectilinear Di stance * In most machine locat?'.on 
problems, travel occurs along a set of aisles arranged in a 
rectangular pattern parallel to the walls of the building# 

In such a situation, the appropriate d5.stance is variously 
referred to as the rectilinear, rectan,f;iilar, or metropolitan 
distance* We choose the fomer and define the rectilinear 
distance between the two machines as 


d * 


- X2 


4 * 




Rectilinear distance is appropriate in some urban 
location analyses where travel occurs along an orthogonal 
set of streets* Additionally, a number of offices employ a 
rectilinear set of aisles and hallways to facilitate the 
travel of personnel* 

(iii) Squared Eucl ide an M stance . In some facility 
location problems cost is not a simple linear function of 
distance. As an example, the cost associated with the response 
of a fire station truck to a fire is expected to be non-linear 
with distance* Depending on the location problem, the cost 
function can take on a number of different formulations, 
xhe most common non-linear form of the cost function is to 
assume it ,to be proportional to the square of the Euclidean 
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distance, or to tafee the distance itself as 

d = (Xj_ - X2)^ + 

There is another reason also for interest in 
squared Euclidean distance. The study of the squared Eucli- 
dean distance problem lays some of the groundwork for the 
Euclidean distance problem "Which is more important but much 
difficult to solve as compared to the foimer. The methods 
of solving Euclidean problems are always of iterative na^ture 
and require some initial approximation. It is customary to 
use the solution of corresponding squared Euclidean problem 
as this initial approxfjaation. In some cases the solution 
of squared Euclidean problem can be used as the solution of 
Euclidean problem also with modera-be error (40), 

D iscrete Spa ce Location Problems - The discrete 
space location problems in^volve the location of new facilities 
at some prespecified locations only. Thus, we have to choose 
the best possible locations for the new facilities from a 
given set of discrete locations. Problems of locating large 
warehouses, factories, plants etc. are some of the tjnDical 
examples belonging to this class. These big facilities cajinot 
be located an^n-Jhere in "tlie plane because of the large floor 
area requirements and/or other fac"tors, 

' L ocation Pr ob lems on Graphs - ILP* s in which we 
can locate new facilities ohiy on the vertices or arcs of some 
prede"beimined netwoidc are included in this category. Problems 
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concerning the deteimination of best locations for facilities 
in a metropolitan city on its road network, in a region on 
its railway neti'Tork or in a telephone sj^^stem on its communi- 
cation netvT'orkj etc* , are some of the typical examples* The 
location problems on graphs are the most practical problems 
because transportation or communication almost always talces 
place along some predetermined network, 

2*2 C lassi fic atio n Based on Ontimization Criterion s 

In this classification we group the location 
problems according to the objective function to be optimized. 
The tvio most commonly used optimization criteria are mini-sum 
and mini-max* Besides these two criteria there are many 
other special criterion which are instituted based on the 
nature of the location problem. For further discussion on 
this classification, location problems have been classified 
as : 

a) Mini- sum location problems, b) Mini-max: location 
problems, and c) special location problems, 

a) Mini-sum Location Problems - In this category 
we place those problems in which cost of total materials 
handling between new and old facilities and between new 
faci?.ities themselves is to bo minimized* This is an 
obvious objective function and most of the location problems 
arising in private sector are solved using th3.s criterion* 
Examples include location of a warehouse for supplying 
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coniino ditle s to custoEiersj location of a factory to supply 
1400 ds ill various markets and location of a new madiine in a 
i'la chine shop, etc, 

"b) M ii ii-»jQa,x.Lp cation P.rphl_em_s ~ In this class of 
prohlems the minimization of the maxiiium of all distances 
(or costs, travel t 3 mes, etc.) from the new facility'' to 
various old facilities is used as the criterion of optimiza- 
tion, The problems of locating a hospital in rural area, 
the placeDiant of firestations in a city, the dotermins-tioii 
of the best location of a police station in a town are 
hand 3 .ed using the mini-max objective. In these cases it is 
not fair to minimize the tota .1 man-miles travelled per unit 
of time, because a minimum man-mile solution may result in 
a situation in which some people, may be very few, have to 
travel a great distance to reach to the service centre and 
when some emergency arises, they cannot be served in time. 

In this case, therefore, we should minirize the maximum 
distance any person is required to travel to reach the nearest 
service centre. Hence mininnax criterion must be used for 
such problems. In general, mini-max criterion is used in the 
location of public sector facilltl.es, specially those which 
are required at emergency. 

Another justifl-cation for minimizing maximum rather 
than -the total distance in the location of public sector 
facilities is that one who installs such facilit.ies does 
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not pay the transportation costs unlUce to the loco-tion of 
plants and warehouses etc. where all materials handling 
costs are to he paid by tho owner of the warehouse. In the 
case of public facilities materials hand.ling costs a.re pa.id 
by those who use these facilities and hence we should mini- 
I'llse the max-inum of the individual costs. 

c) ~ 111 addition to the 

tifo optjmization criteria dj.scussed above ^ there are some 
location problems in which we are forced to optimize using 
some' need based criterion. For example, we sometimes come 
across the problems of netifork design in which several points 
on plane are given, and we are required to join these po.iiits 
with smallest possible total li;t; length, These points can 
be joined in two ways. One .i.s without creating a.ny neir 
nod.e which means that junctions are only at given points, 
a,nd the other is when this restriction is released, i.e,, 
we can create new junctions also. When the second method, 
which i.s more realj.stic, is used we have a lo ca- tie n problem 
in which We will have to deteiinine tho optimum number of 
extra junction points, their locations and then the optimum 
l.lnlring betij’een various points so that the total linlc length 
is minimized. Those problems arise in real life when we 
have to join different towns by minim.um possible transpor- 
tation nettrork or different points in a X'fater supply system 
by mindmuia length of pipe line etc. 
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Anotlier spec^.al location prollom arises while 
deteiiiiining the locations of einergoncy service facilities* 
For such problems it is desired to locate the facilities 
such that no requirement centre is farther aifay than a 
prespeeified critical distance or critical travel t.Uiie from 
at least one facility, i*e,, all the requirement centres 
are “covered”* Thus we want to determine the minimum number 
of facilities and their locations so that the objective is 
accomplished (total covering problem), or given tha.t a 
particula.r number 03!* facilities are to be located, to deter- 
mine their locat.ions such that ma^cimum number of require- 
ment centres can be covered (partial covering problem), 

(11a,ss i tica.t ip n B a se d on th e llu^j,r a_ncl_ Chcm^ctejr^ 
of the, Facilities to. be 1 

Under this classification, the location problems 
can be grouped into the following c.lasses : 

a) Single fac3J.it3'' location problems, in which 
onJ-y one new facility is 'be be located. This problem when 
formulated in continuous space with a. mini- sum objective is 
cal 3 .ed the "Depot Location Problem", The depot loca,tion 
problem with Euclidean distances is called the "Steiner- 
Weber Problem**, *’ileber Problem", or the "General Feima.t 
Problem", 

b) Multi- facility location problems, in which we 
consider the location of ttfo or more new facilit3.es- 3 ssent 3 .a 



there are t^-xo types of multi- facility location problems : 

(i) Problems of locating dissimilar facilities j 
i*e, those facilities vhich serve different purposes (called 
DSPLP* s henceforth)* 

(ii) Problems of locating similar facilities, i*e, , 
those facilities vrhich serve the same purpose (called SILP* s 
henceforth)* In SILP^ s we have to deteimine optimal alloca- 
tj.on also, i,e,, the allocation of new facilities to serve 
the various demand centres. 

DSILP* s can again be divided according to how tiie | 

old and new facilities are interactiong, i.e,, vrhat is the \ 

pattern of materials handling betvreen them. Using this f 

criterion, we can subdivide these problems into three categoriess i 

A* When new facilities interact with old ones only. | 

B, When new facilities' interact within themselves onl7/. s 


C, When new facilities interact with old ones as 
well as within themselves. 

It seems, at first sight, that the new facilities 
in any DSILP must ha.ve interaction with the old facU.ities 
as well as within themselves, i.e., they must belong to the 
third category only because if there is no interaction betvreen 
the new facilities then a multifacility location problem in 
which n new facilities are to be loca.ted can be converted 
into n single facility locati-on problems. On the other hand, 
if there is no interaction beti/een the new and old facilities 
then there is no location problem (all the new facilities 
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must be located at the same place and this place nay he 
any\rhere in the plane 1)* Hence only the prohlems of the 
third category seem to he meaningful and this classification 
seems to he needless. But when we are locating the facili- 
ties in discrete space or on networlcsj then almost alimys we 
have a restrj.ction that no more than one facility can he 
located at one given site. Under such circumstances, an 
indirect type of Interaction develops between the new facili- 
ties, This validates the above clp.ssification. Of cotirse, 
location prohlems in continuous space need not he classified 
in this marine r. They alwa^/s belong to the third ca.tegory 
only. In discrete space, the proh3.ems of first and second 
category are referred -feo as the ^Assignment Problem" jand 
"Assigranent of Facilities to Locations (or the Quadratic 
Assignment) Problem" respectively. The problem of third 
category j.s yet to he posed and. solved* 

SPLP * s can he further divided 5n t\>ro ways. Following 
ti/o t3?pes of prohlems emerge out from tho first of them : 

A. Prohlems in which the number of new facilities 
to he located is known. Such prohlems arise in the location 
of. public sector facilities where the number of new facili- 
ties is limited due to financial constraints, 

B, Prohlems in which the number of new facilities 
to he located is unknown, * These are prohlems arising in 
pri.vate sector where the optimal number of new facilities 
is to he determined in addition to their locations. It 
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seemsj at first sight, that in this case also we shall have 
restriction on the nuiiber of neiT faci3.ities due to financial 
limitations r''i hence it must he known. But here the diffe*-- 
rence is that if we decide to locate too less riumiDer of new 
facilities then the saving due to this wj.ll reflect in the 
increased transportation costs since all the costs are 
charged to the owner of the facjJLities. Hence in such 
situations we have to determine the optimum numher of new 
facilities such tha.t the overall cost of instal3.ing the new 
3 ?acilitios and tra.nsportation is minimum. This proolem, 
wlien formulated in continuous space with Euclidean distances, 
is called ’'Location-Alloca.tion Problem'* and in discrete 
space, the "Plant Location Problem" or the "Warehouse Location 
Problem”, 

The other \ra 3 r of dividing SELP^s gives following 
ti -70 tyjoes of problems t 

A, Uncapacitated SI^LP* s, in which there is no 
upper limit on the capacity of the new fac 5 J.ities to be 
lo cate d* 

B, Capacitated STLP* s, in which there is a restrict* 
tion on the capacity of tlie new facilities, so that not more 
than a given amount can be supplied from any new faciJLity. 

Ifo capacitated SILP was formu3-a ted in continuous 
space till very recently* In 1972, Cooper (7) formulated 
and solved the problem and called it *^he Ti’ansportation- 
Lo cation Problem", which actually is the capacitated SELF 



iii contjjiaous space or the ”Capacitated Location- Allocation 
Prob.lem", 

SELP*s helonging to the above tivo types In discrete 
space are called 'incapacitated Plant (or ¥a.re house) Location 
Problems'' and "Capacitated Pla.nt (or Ifarehouse) Location 
Problems" respectively, irhen the optimal number of new 
facilities also is to be determined* 

Fig* 2,1 depicts the classification on the basis of 
number and characteristics of new facilities in a concise 
format, 

the. PJ-aryning. Hpi^zoR * 

tfe can also classrify the location problems as to 
whether it is 

a) a static location problem 
or b) a dynamic location problem. 

In the first category we keep those problems in 

which the input data, i,e, the parameters of the problem (the 

numher and location of demand points, their level of demands, 
etc.) bolong to only one point in time, Thl.s implies that 

the planning is done as if the parameters of the problem are 

not variable with tiime. The second category includes those 
pioblems for which parameters are functions of time vrithin 
the planning period. 

Every effort has been made to envelope the whole 


range of the literature presently available on modem loca- 

■i - 

tion theory in jthe proposed classifJ.cation scheme. However, 
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Fig, 2,1 j Classification of Facility Location Problems on 
the basis of number and characteristics of new 


facilities 











n o 
-Lo 

tho literature on the location problems is expanding at a 
very fast pace and it is becoming increasingly difficult to 
keep a complete tract on them. Therefore, there may be some 
location problems which have not been included in any of the 
categories suggested above. It is not to say, therefore, that 
these classifications are complete in themselves. Still 
the3r give a broad idea of the entire domain of modern location 
tliGory and help to overview it. 



Chapter 3 


SIH( g.E FACILI TY LOCATION P BDBL MS 

The objective in single facility location studies 
is to detemine, for the new facility} an opt5nal location 
which will minSnize appropriately defined total materials- 
handling cost or total distance function. The cost of 
materials handling is considered to he proportional to the 
distance. 

There are a numloer of real 3.ife situations which 
involve the loca.tlon of a single ne%/ facility in an existing 
layout. Some typical examples of single facility location 
proh].ems are the location of a new lathe in a manufa.ctur3ng 
job shop; a new warehouse relative to production facilities 
and customers; hospital, fire station, police station, or 
library in a metropolitan area or a new power generating 
plant etc, 

3,1 Single Facility Lo cati o n J'todel^ s : 

In literature, the problem of single facility 
location is posed as follows s n existing facilities are 
located at Icnown distinct points a new facility 

is to be located at a point X; costs of transportation are 
incurred that are directly proportional to an appropriately 
determined distance between the new facility and existing 
faclJLities, Let d(X, P^) represent tlie distance travelled 
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per trip TDetv-/een points X and P^, If represent the 
numher of trips made per year between the new facility and 
existing facility i, the total annual cost due to travel 
I' betvreen the new fac5.1ity and all existing facilities is 

given by 

m 

f(X) = I wj d(X, Pj_) (3.1) 

i=l 

The w^ terms are referred to as "weights”. The 

r 

single facility location problem is to determine the loca- 
tion of the new facility that lafnimizes fCX), the total 
annual transportation cost. 

' Depending upon the type of distance function 

d(X, Pj_) used in a particular problem, the single facility 

location problem is classified as : 

1 , Eectilinear distance location problem, 

2 , Squared Euclidean distance location problem, 
and 3. Buclidean distance location problem. 

The formulations and solution procedures for the ^ 
problems are presented in the following sections, 

SaIj 1 ^c t-ilinear D ista nce locati on Pr oblem — The recti- 

linear distance location problem can be stated mathemati- 
cally as given below. If the co-ordinates of the new 
facility are (x, y).and of the existing facility i are 
(ai, b^) so that X = (x, y) and P^^ = ^ 1)5 location 

problem 'is given .by v'," 
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ra 


Min 

x,y 


f(x, y) = I (Ix-aj^l + jy-bij) 
1=1 


(3.2) 


The problem given by (3.2) was first solved by 

Francis (10) in 1963, He gave a stople procedure to find 

the optimum values of x and y co-ordinates of the new 

facility and discussed some interesting properties of the 

problem. He has also given an algorithmic procedure to 

1 

construct the contour lines for the problem. His solution 
procedure is based on the duality theory of linear progra- 
mming. 

3,1, 2 Squa re d ^c lj dean D l^stance, 1 p^a tjp n^ Pj^b^lej^ - The 
single facility squared Euclidean distance location problem 

can be stated as follows s 

m 


Min 


f(x,y) = 1 


(x-a^)^ + (yvbi)^l (3.3) 


i=l 

An interest in the squared Euclidean dj-stance 
problem or the gravity problem (another name of this 
problem whose justification will be given in a minute) is 
no t new , Acco rding to Kelly ( 20) , inte rest in thi s p roblem 
dates back atleast to the work of Lagrange, McHose (29) 
has used the basic concepts of dj-fferential calculus to 
obtain the optimal solution for the gravity problem. 


A contour line is a line of constant cost in the plane. 
Locating the new facility at any point on a contour 
line results in the same to^l cost. 
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m 


21 


f(x, y) = I (jx-a^l + jy-bij) (3.2) 

1=1 


The problem given by (3,2) was first solved by 
Franc.is (10) in 1S63, He gave a simple procedure to find 
the optimum values of x and y co-ordinates of the new 
facility and discussed some interesting propertl.es of the 
problem. He has also given an algorithmic procedure to 
construct the contour lines^ for the problem. His solution 
procedure is based on the duality theory of linear progra- 
mming. 


3 ,1, 2 Squa re d ^c lideari Dis tance Lp cation PTOb^-em^ - The 
single facility squared Buclidea.n distance location problem 

can be stated as follows : 

m 


Min 

x,y 


f(x,y) = ^ ^ w^[ (x-a^)^ + (yvbj^)^] 


(3.3) 


An interest in the squared Suclidean distance 
problem or the gravity problem (another name of this 
problem whose justification will be given in a minute) is 
not new. According to Kelly (20) , interest in this problem 
dates back atleast to the work of Lagrange, McHose (29) 
has usod the basic concepts of dj.fferentioJL calculus 
obtain the optimal so3.ution for the gravity problem, 

A contour line is a line of constant cost in the plane. 
Locating the new facility at any point on a contour 
line results in the same total cost. 


.1 
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Recently Cooper (4,5,6), Knhn and Kuenne (21), Mlohle (30) , 
Weiszfeld (41), Vorgin and Rogers (40) have used the solu- 
tion to the s.ingle facility gravity prohlen to initiate a 
search for the optimum solution to the single facility 
Euclidean distance prohleni or the generalized Fermat problein(21)o 
The gravity problem has a simple solution. The 
optinal co-ordinates of the new facility are the weightad 
averages of the corresponding co-ordinates of the existing 
facilities or in other words the centre of gravity'* so3.ution 
is optimal to this problem;, Hence, the title gravity problem. 
Contour lines for tliis problem are concentric 
circles centered on the optimum loca.tion of the new 
facility, i,e. , the centroid of the exisWJig facilities* 
locations, 

3.1,3 Bu cli dean Distance Location Problem. The Euclidean 

problem may be stated as : 

m 1/2 

r(x,y) = I Wit (x-ai)® + ] (3.4) 

i=l 

A version of the IJuclidean pTOblcsii for the case 
m = 3, Wj_ s 1, for i = 1,2,3, was posed, purely as a 
problem in geometry, by Feimat early in the seventeenth 
century, a,nd was solved by Toricelli prior to 1640. The 
problem was studied by Steiner, a Swiss mathenatic.ian, in 
the nineteenth century, and byl'Jeber, a Geiman economist, 
early in the twentieth century. Interestingly enough. 
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however, it was not until tlio %-rorlc hy Kuhn (21), in 1962, 
that the problem could be considered to be essentially 
completely solved* 

The -Tuclidean problem can either be so3.ved by the 
modified gradient iterative procedure of Kuhn (23) which 
is guaranteed to converge to the optimal location; a 
discussion of convergence for this is given by Katz (19), 
Kuhn (22), and He is zfe Id (41); or by the Hyperboloid 
Approximation Procedure (HAP) given by Syster, Hhite, and 
Hieri'/ille (8)* HAP is an iterative procedure based on the 
use of an approximating function and can be used to solve 
the rectilinear distance location problem as well. 

Unfortunately, exact methods for constructing 
contour lines are not avaiJ.able for the I&iclidean problem, 
except for the simplest cases where there are only one or 
tv7o existing facilities* For other cases, a relatively 
simple method is to obtain approxiiiia,te contour lines by 
evaluating the cost function over, say, a rectangular grid 
of points covering the ranges of (x,y) values of interest* 
The contour lines can then be drawn by interpolating 
between grid points* Alternatively, one ca.n assign a 
given value k to fCxjy) in (3,4), pick a value of x, and 
search over y for the two values that yield the value k* 
The process is continued for suceessi'^' values of x until 
a, family of points is obtained fcr the contour line having 
value k. By repeating the same procedure for different 
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values of Ic many contour lines ca,n be drawn* 

, Fj’PPospA JL?jte n§3p n, pX t]*®, • 

In certa.iii classes of location problems, the 
munber of existing facilities or demand points may be exce- 
ssively large, ThJ.s would cause high computational effort 
for the optlna.! location of the new facility to obtain. 
Further the rep3:esentation of la.rge number of demand 
points or facilities becomes inconvenient, , For example, 
if a p 3 ?oblem involves the location of a new facility 
(distribution point) in a densely populated urban area, 
it would be most difficult to represent the demand points 
(each resident) with a point, 

A common method for solving such problems is to 
divide the total area into sub-areas and to represent the 
population of each sub-area as an aggregate **point**, i*e*, 
the centre of gravity of that sub-area, A shortcoming of 
this method is the error introduced by the assumption 
that the distance travelled from any point in a sub-area to 
the new facility equals the distance travelled frcsa the 
centre of that sub-area, to the new facility, models 

which are more accurate representa.tion of -the actual 
situation are proposed here. 

In the first model, we divide the locations into 
groups, each group having seme average population density, ■ 
The second model utilizes a demand density function defined 
all over the region under consideration. We shall 


25 


describe the details of these models one by one. 

3.3 Model I .Rect angg il a .r Are a An n roach : 

In this model the total region under considerat5.on 
is divided into rectangularly shaped areas of known dimen- 
sions* Each area is assumed to have a uniformly distributed 
population density, i.e*, the demand density at each 
point in a particular area is same. The dememd density 
maybe different for different rectangular areas. ITaturally, 
the accuracy of representation inc rea.se s with higher nuinber 
of rectangles used for dividing the total region* 

The idea of using uniformly d.3.stributed populations 
for ind.i vidua,! rectangles is not a new one. Several contri- 
butions to the facilities location literature that use 
this concept have appeared. Learner (25) used heuristic 
methods and straight line distances to allocate sources to 
the interior of a,n area with uniformly distributed densities, 
Francis (13,16) has described several models with contljiuous 
destine, tiohs. Bender and Goldman (3.) and Witegall (47) 
have published models of post office systems where the 
districts have unifomly distributed popu.lations. Hewman(32) 
has proposed a parking lot location problem and has given 
a elution for the special case where the facilities are 
straight line segments* A solution method for the case 
where distances are given by the generalized Ip 
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1 

function lias not •'/et been developed (36), 

In order to fomulate the total travel cost, 
consider a t^dPical rectangle with length b^ - and 
width d^ - Cj, Let be the location of the new facility 
in the region under investigation. If represents the 
deriSit 3 r for the ith rectangle, then the travel cost for 


this rectangle is given by 
^i '^i 

(X^^) =111^ / / dCX^^, Z) dy (Ibc 


“1 






where X = (x, y) is any point within the rectangle. 

Let n be the total nunber of rectangles. Then the 

total travel cost for a.11 the rectangles will be 

n b. " 


n 


i 


f(Zn) = I f^U-y) ^ I m. / / dCXn ,X) dy dx (3.5) 

1=:1 1=1 

which is to be minimized with respect to the location of 
the new facility, i.e,, X^ (x]_, y^)* 

We shall solve this problem for all the three types 
of distances discussed ea.rlier. 


^ If we are given two points in plane, Xp == and 

Xg =* (Xg, y 2 ), then tte 3^ distance between Z 2 _ and X 2 
is glvan-by , p ^ 

=t| xi=s:2| + |5l-y2| 3 

Bectllinear and Euclidean distances are two special cases 
of Ip distance with p=l and 2 respectively. 
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Dj.stanc e j;. o cation _P rolDlem - When 
rectilinear distances are used, the location problem 


he cones 


n 


Min^ fCx^jYl) * I / / tl x-x^j+ jy-yj] dy cbc (3,6) 

1^1 8 , « C 4 

i=l ^ ^ 

Fortunately, (3,6) is separable and the optimum 
values of and y^ can be found independently. This follows 


DO cause 


Td^ di 


n b , d. 

1 X 


- I “1^ J J Is: - x^l dydx + I f f }y-yj djt 

1=1 a, ®3. 


I (dj-op / |x-Xj^|ax+ I nj^Cbj^-ap / |7- “ 

is=l " ^-i i=l 


J % / |s:-Xj T f |y-yj dy 

i=l i— 1 c. 


= fiCxj,) + fgCyi) 


(say) 


where = m^^ (d^-c^) and - m^i^ Cbj[-a^), 

The optimum x^ and y^ values can therefore be 
found out by considering the location of facilities among 
population lines along x and y axes rep sect! vely, iJe shall 
discuss tl^' precedhre for miii^izing fiCxj^) only because 


tlie mlnojnization of ^ 2 ^^^ done by following the same 

procedure. 

The minlaization of f^(x 2 _) can he thought of as a 
one dimensional location problem,. In this problem we have 
to determine the opt.lma,! location of a single facility with 
respect to line destinations placed on an axis. These Ijiie 
destinations can overlap each other. It is convenient to 
pose an equ3- valent problem whore none of the destinations 
overlap. For example, if there are two destinations 
[a^ bj_] and[a.j, b^] , such that a^ < a^ , S’^d bj > b^., 

then 



and three non-overlapp?.ng lines [a^,aj3 , ta-j, b^l, and 
[ b - , b. ] are obta.ined. After the overlaps have been iden- 
tified, there will be n* line destinations where n* ^ n, 
¥e can arrange tie destinations along tl^ axis from left 
to right and label them [ r^. , s^], i=l,.,, ,n^ where 
r^+1 > s^ , r . < s - , at - r» , b! = s . and the constant is 
A^ , ¥ith these modifica.tions our objective function w5.11 


look something l.lke 

n* 
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All the terns in convex 

is also convex, A teim such 3 -s A^ 


, and hence 

/I |x-x^|ax 



derivative of that term* 


From the figure it is evident that %fhen x, < al . the 

' ' J-'' ' 1^'' ' 

slope of the term corresponding to £r^, s^J is a negative 
constant and vrhen > h^ tl^ slope is positive value of 
that same constant* It i|^ PQss^|>le 'to evaluate, ! f^(x2^)/dx3_ 
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I . . 

at the values j i = 

Let the absolute va.lue of the constant slope 

corresponding to [r^j s^J be t^, i = , Tlion 

t-. = At (bt-at)« Whence-, < r., 

2 . 1 1 11 

d fiC 2 C 3 _)/<ix 3 _ = ^ I t^ = - M (say) 


i=l 


Consequently 


d f^(x2_)/dx-j_ 


S^ = - 2 -I t^. 

D=1 


(3.7) 


Since f^(x^) is convex, this expression enables us to find 
the region in which f;L(x^) is a miniriium. For this purpose 
(3.7) is evaluated for successively larger intervals i, 
until it becomes either zero or positive. If it becomes 

f ^ I 

zero for some value of i, say i* , then Sn_ ± ± ^i+1’ 

where x^ is the optimum value of x-j_. If (3.7) becomes 
positive for the first tjjae when i = i* , then r^ £ x^Als* . 

In the la.tter case , the exact pos.ltion of x^ can be found 
o’ut by using the derivr.tive plotted in Fig. 3,1, whereas 
in the foimier case x^ can be taljen any wliere in the .Interval, 

}|5 

The co-ordinate y^ of the op t3mum location of the 
new facility can be found in a similar manner* 

We have previous3-y discussed tlie importance of 
contour l.ines while deal.ing with the discrete demand 
fac3J.ity location problems. The contour lines have the 
same importance In continuous ctoand locat3.on problems also 
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and no\T method 'Will be described to construct them for 
continuous domanid location problem. 

In this case the equa.tion of the contour line will 
be 


k 


n 

h. 

^-i 

Z ">1 

/ 

/ 

i=l 

a. 

c. 

1 


where k is a constant. 

The above equation can be rewritten as 


n 

^i 


n 


% 

r* 

1=1 

(d.^— C jj) J 1 x- 

" ■■ Si 

-X3_ Ite + 

i 

i=l 

m^_(b^-a^.) 

/ 

C. 

1 

n 


n 

"i 



I 

J l3:-x^| dz + 

Z h 

/ 

! y-yil <ay 


i=l 


1 

H 

iii 

®i 




a 


How, suppose that the point (x 3 _,yq_) is such that 
® Tlien we sha.ll have 


n ^1 ^ 


yi 


k = I I (x-x^)doc] + I / Cy^-y)l 

. i=l 


a=l ^ C-- 


^i 


+ / (y-yn)dy] 


Rear ranging the above equation, we obtain 


o o 
ad 


n 


n 


n 


k 


= ( I +( I B )y^ - [ I (a,.+l) )A. 

i=l i=l " 


i)A 3 - 


n 


XI 


2 . 2. 


[I ( c^._+dj_)B^ 3 jj+-l I [ Aj_( a^+b^)+Bj ( + dp 3 


i=l 


i=l 


This is a second degree equation heti’reen and 
and is the required contour line for the rectilinear 
prolDlem* This equation is valid for the rectangle -which 
extends to [a^, h ^3 in x-direction and to [ Cj^jd^. 3 in 
-direction. One such equation is to he deteiTAlned for 
each of the rectangles to construct one contour line. By 
taking different values of k, many contour lines may he 
constructed. 


_3..3. 2 Squared, guclidean J)lstejice_ Pjph lein ~ The 
continuous demand squared Buclidean distance location prohle 
or tho continuous demand gravity prohlem - 1/111 he given hy 


Min 




f(xi,7i)- I n 3 _ / / [Cx^x )^+ Cy-y-,)^}3.y dx 

i=l a. 

.{3,8) 


®-i ^i 


It is oh served that ( 3,8) is separable in x^ and y^ and 
we can minimize it hy breaking it into -&/o components, as 
shown helow^ i ^ 

The component corresponding to x. 5-S 

” ^i i 


n 


h 


.! “l / / dy dx 


i=l 


^i <^1 



n 


0-0 

OO 


’’i 

y f 2 

^ els 

i=l a._ 


Differenti3.tion %-J'ith respect to s^ gives 

2 

^ 1 /flTT-. = y A / r •y^'v- 1 
1 ' 


n 



I 

^i 

/ 

i=l 


^1 

11 


^i 

1 

A. 

/ 

i=l 


a. 

1 


For obta.injiig a niinj-mum V3.1ue of f^(s 2 _)} "'f^e uust equate 
df^(s^)/dx^ to r.ero, getting 


n 

I i A. ;/ -2 (x-St ) dx = 0 


1=1 


a 


1 


Solvi.ng the above equation for x^, vie have 

n n 

2 2 

~ ^ h 3 / [2 I A (b^-a^) ] 

1=1 


1=1 


Similarly optninim value of y co-ordinate can be determined. 


Go ntour L lnes» In this case the equation of contour line 
becomes 


n 


n 


k = I A. / Cx-X 3 _)^d 2 c + " I J (y**y^) dy 

b i=l ^ ^ i-1 ■ 
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mvaluating the integrals and rearranging, ve ohta.in 
n n 

k =[ I + [ I B^(a^-c^) ] - 


n n 

[ I I - i I ®i^'^i~°i^ 3 + 

i=l 


+ I I [A^(h^-a^) + B^(d^ - c|) ] 


which is the required equation of the contour line. 


3 ,3. 3 Buclid e an p i t tance. Lp.cation Problem - In this case 


the problem will he given by 

n b, 

Min ■ r „ f 


\ % 


^ y f(x3_,y3_) = I / / [ (--3^) + (y-yi) ] ^yc^:x 

“■1’ 1 -1=1 ' a. c. 

^ ^ ^ (3,9) 


For a minimum value of f(3r2.}yi)} '''^© must equate 
the partial derivatives of f^(x;j_,y^) with respect to 
and y-^ to zero, obtaining 




n 

^i 

% 

9f/ 

ayx = 

I “i 

f 

/ 



3=1 

^i 


and 
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% 

3 f/ 

11 

1 m, 

1=1 ^ 

/ 

H 

7 

°i 


3 _ /c^ [(x-x^Li^+Cy-y^)^' 


dyd2?=0 (3,10) 


(y-yx) 

_ a2./. 


,2 1 2 


dyax=0 (3.11) 


JL 

p 2 ^ 

Let [ (x-X 3 _) + (y-y^) ] = d, then (3,10) and 

(3,11) after rearrangement give the following iterative 


formulae % 

n - 

I \ ! J [z/d^^^] dydjc 

(hi-l) i=l °i 

^1 " h^'X' 

I / f [1/d^^^] dydx 

i=l a. c. 

H- II. 


and 


n 


■b^ 

1 1 


(b) 


I Ei^ / "/ [ y/d ] dych 


y-i 


(h+l). i=l 


a. G. 

1 1 


n 


I r 

is=l 


y 


[l/d^^^^] dydx 


a. c. 

1 1 


(3. 


The superscripts in the above formulae indicate 
the iteration number. 

The initial values of and y^ in (3,12) and 
(3,13) can be tafeen as those optimal for a corresponding 
squared Suclide an problem. The stopp5*.ng criterion may be 
based on successive changes in the x and y co-ordinates of 
the new facility or on successive change in the value of 
objective function or on both, Sinc^ our objective func- 
tion is convex the iterative procedure is convergent. 

The minimization of (3,9) can be done by nonr- 
linear programming also and it is interestjng to compare 
(3,12) and (3,13) with the iterative procedure used when 


this prohlGia is solved as an unconstrairBd non-linear 
programming prohlem using the steepest descent method. 


From (3,10) we get 
^ (h) 

^ ^i / J ^ ^ 

i=l 8 .^ 


n 


"^i “^i 


I “i / / 

4 _i a. c. 

3-=l 1 1 


3f/ 3 


(h) 


(3^14) 


Siibstj.tuting (3,14) in (3,3.2) and solving for x- j we obtajii 

n -bj^ _ai 


(h+1) 


"1 




n 


^ To.^ ^ ^ [l/d^^^ Hjdx: - 3 f/ 

i=l - - - ^ 


^i *^1 


I 7 ^ / 1/d^^^ ] Uydx 


ir=l 


w 


(h+1) (h) Xh) " ^ (h) 

or - t 9f/ 

n ^ 

where J m^^ / / [ 1/d^^^] dy cbc } 

a ^4 
i=l it 

Similarly 


■1 


(3,15) 


(3,1S) 


(h)_^(h) 


(3.17) 


It is evident that our proce dare differs from the 
steepist d® 2 cent method since the value of t^ gj.ven by 

(3.16) need not be an optimal step length. However, the 

{ '■ 

ease with which the value of u ' is calculated using (3,16) 
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results in a solution procedure uliicli is easier than the 
steepest descent raethod. Furthermore, no 'background in 
non-lihear programming is required, to apply our procedure. 

Contour Lines , S.lmilar to Suclidean distance dj.screte 
facility location problem, .in this case also no direct 
method is available to construct the contour lines. They 
are to be drawn as discussed previously for the discrete 
djemand problem, 

3,4 I-iodel II .:....D ensi.ty:. .Funct ion Appr oach : 

The second model to be described defines a 
demand dens.ity function g( 2 :,y) all over the regixm under 
consideration. This means that the demand at an infixii.te- 
simal area djxdy at point (3c,y) is g(3c,y) d:!cdy, ¥e shall 
solve the single facility .location problem using this 
approach for all the three distances prev.iously discussed. 

It will be assumed tha.t the region und^er consideration 
extends to a,b in x- direction and to c,d in y~ direct ion 

3,4,1 Bectili n ear Dist ance Locat.ion Prob lem ~ The total 
cost of transportation in this case, ¥.111 be given by 
b 

^‘(2Ci,yi) = ^ ^ g(2:,y) [ jx-x^l + jy-y-jJdydx 

b ^d ; ■ b d 

= / / g(x,y) dydx + / / gCxjy) j y-yj dydx 

.a ' ' ,''0 " ■ , , ■ a O' y' ‘ , 
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results in a solution procedure which is easier than the 
steepest descent method. Furthermore, no hackgrouiid in 
no'n-lihear programming is required to apply our procedure. 

Contour Lines * Similar to Euc3-idean distance dj.screte 
facility location problem, in this case also no direct 
method is available to construct the contour lines. They 
are- to be drawn as discussed previously for the discrete 
demand problem, 

3.4 Model II D ensity fian ctio n . Appr oach t 

The second model to be described defines a 
demand density function g(x,y) all over the region under 
consideration. This means that the demaxid at an infi.xii.te- 
simal area cbcdy at point (x,y) is g(x,j^)d3cdy, ¥e shall 
solve the single facility location problem using this 
approach for all the three distances previously discussed. 

It will be assumed that the region under consideration 
extends to a,b in x-direction and to c,d in y-direction 

3,4,1 Rectil i near Dist ance Location Prob lem » The total 
cost of transportation in this case, will be given by 

. - ^ - ::;h : b ,4 : ' : _ 

f(xi,yx) = ^ ^ gU,y)[ jx- 2 ^| + P dydx 

a c . 

b d b d 

~ j f g<3C,y) I J g<3:,y) iy*-yjdydx 

a 'C , . . 1 ' ■ a ^ e' - 



t) d d b 

/ Ix-X^l J g(x, 3 r) dydx + / 1 y-y^| / g(x 5 y)dxd: 


a 


d 


= I I F(x) cbr + / |y-yj G(y) dy 

a c 


where F(x) and G(y) are the ma.rg3.nal dens3.ty flinptions of 
the den.s3.ty faiiction g(x,y) with respect to x and y 
respectively, 

S3.nce the two terms in the express3.on for f(x,y) 
are .independent of each other can minjmize them separately. 
Let fCx^^jyi) = 

>1 h 

then f^( 2 ^) I (3:2_-x) FCx) dx + / (x-x^) F(x) dx 

a- ^ 

= / x_ F(x) dx - / xF(x) dx + / F(x)dX“ 

a a \ 

h 

/ X_ F( x) dx 

Differentiating with respect to Xx g®t 

h :■ ^ ^ 

d f^(xj_)/dx^ = X^Kx^)+x^F(xj_)+ / F(x}dx- / F(x)dx t 

xx a 


ic^l -F(x^)-P(x^) i 
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13 

^ F(x)(i;c 




^ F(x) dx 
a 


For value of f^(x^), this must he equal 

to zero, so that \ro get the conditil.on of optimality as 


1 ^ ' 

J fix) = 

a 

h 


h 


j F( x) dx 


Sirj-ce j F(x) dx = 1, hotli the two sides of the above 

a ■ , 

equation must he equal to 1/2, i,e* 


/ F(x) dx = 


h 

/ F(x) dx » \ 


The point x^, satisfying the above property, is 
called a median point. Hence we conclude that x-coordinate 
of the optimal location of the new facility is at the 
median of the marginal density functden with respect to 
X, Similar result applies to y- coordinate also, 

How it reman.ns to detem.lns the median of a given 
probability density' flinction, say f(x). This can be done 
easily as described below. 

Suppose f(x) Is defined over [x^, xj , i.e,, 

/ f(x) dx 1 f 



Also let m te the median of f(x). Then 
m ' ^2 

/ f(x) dx = / f(x) dx = -i- (3,18) 

^1 ^ 

Let / f(x) dx = F(x) 

Then (3,18) can he written as 

F(n) - F(x^) = F(x 2 )~F(m) = i [F(x 2 ) -* F(x^)] 

or 2F(m) - F(x^) ^ ^(xg) = 0 

This is an algebraic equation in m which can he 
solved, at least hy a numerical procedure, to g:lve the 
required median of the probability density function f(x), 

Coir'coiir lines . The equation of a contour line i/itii total 
cost Ic will be 
b a 

!: = / / g(x,y) [ |x-x^| + ly-y 3 _|] dy dx 

a c 

b ^ 

= / [x-xj FCx) dx + / ly-y^^l G(y) dy 

(x»-x^)F(x)dx + 

d ' ^ 

/ Cy-y^)G(y)dy 
% ■ 


a 

a 


c 

b 


(xt-x) F(x) dx + / 


/ (y -y)cKy)dy + 
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Siiuplifyilig and us3.ng 
b 

/ F(x) (fe = / G(y) dy = 1, 

a ^ 

b 

/ 3cF(x) dz - = raeaii of the marginal density function FCx), 

a 

4 

and I yG(y)dy = m = mean of the marginal density function GCy), 

V 

C 

ire obtain the equation of the contour line as 

b d b 

k = 2 / xFCx)dx +2 / yG(y)d3r+ ~ / F(x)dx + 



y^ - 2y / G(y) dy - m ~ m 

4 "V ^ 


The above equation can be written explicitly in 
teims of and y^ if we know the marginal density functions 
F(x) and GCy), which in turn can be determimd from the 
e::pression for the demand density function g(x,y). 

3« 4,2 S quare d Eucli dean Mstance Location Pr o blem ~ In th3.s 
case we have the following escpression for the total trans- 
portation cost, 

b d 

2 2 

f(xi,y^) = / / gCkjy) { (x-x^).+(y-y^) 3dy dx (3»19) 

a c 

'For getting 'a mihimuBi value ©f f(x-j_)yiL), we Enist 

'if - * ■ 

equate its partial -derivat.ives with respect to nnd y^ to zero* 
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= 0 gives us 

Id d 

/ / gCxjy). 2 (x 3 _“x) dy d^r = 0 

a c 


Solving for x^, we oltain 
"b d 

/ / "x g(x,y) dy djc 
_ s. ' c 

1D d 

^ ^ g(Xjy) dy dx 

a c 

Td 

= / X Kx) dx, since denominator is equal to 1, 

a 


== m 


X 


Similarly y^ = m^ 

Co ntour Lines . The equation of a contour line with total 
transportation cost k vrill be 


1 , „ / -/ g(x,y) [ (x-x ) + Cy-y^) ] dydx 


h d 

“ / ■/ g(x,y) (x^ + x^ - 2xX|_ + y^ + y^ « 2yy^) dy da: 

a G ' 

= n “xa - 2^ “x ■*• 5^ * “Va - “y 

where and m , are respectively the second moments of 
F(x) and G(y)# 
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Hence, ve conclude that the contour lines for 
this problem are concentric circles with centre (m ) and 

2 2 i y 

radius (ni^ + m « - m„ + k) • 

y X2 y2 


3,4,3. Euclidean Distance Location Problem. When Eucli> 
dean distances are used, the e^cpression for total cost 
will be 


b d 


r 2 

f(xi,y^) - I ! g(x,y) [ (x-x^) + 
a 0 

For an optimal value of x-j^ we must ha.ve 



2 


JL 

2 


dydx 


af^(x^,y^)/ 3 x^ = 0 


Differentiating and solvliig for x^, we obtain 
b d xg(xjy) 

j ; — 


"1 


2 2 ^ 

a c [ 1 

b d 

J / dydx 

a c 


g(3c,y) 


[ ( x-X 3 ^) +( y-y ) ] 


2, I- 


Hence wiH have to be found out iteratively, 
i,e, , first we give approximate values to x^ and y^ and 
then improve upon the value using the above ecjuatlon. 
Similar iterative procedure is to be u^d for dote imining 
optimal valuer of y^^* The stopping criterion and starting 
point similar to the previous appioach may be \ised here 
also* . . 
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The use of noixLinea.r prograiiuiiing is also possible 
for solving the Euclidean distance location problem, but 
it uill not be as efficient as the above itera-tr'.ve procedure. 
Contour lines, once again, cajinot be dravrn 
directly. An ind3.rect approach discussed provn.ously is to 
bo used. 


Chapter 4 


IdXTI-FAC ILITy LOCATION PBDg .lS'iS 


In our previous discussion of facility location 
problems we treated the case of a single rew facility to 
he .loca.ted relative to an area having continuous demand for 
the facility. In this chapter, the previous analysis is 
extended to include the prohlera of locating multiple new 
facilities with respect to an area which has continuous 
deraands for each of the new faciliti.es. Thus the single 
facility locat.ion problem can be considered to be a spec.ial 
case of the multi- facility location problem treated here. 


4 .1 . Discrete Demand Multi- Faci li ty L ocatio ii PTObl em s 


First, we shall give a brief account of discrete 
multi- facility location problem and then switch over to 
its continuous version. The discrete multi- facility location 
problem can be foimulated as follows. Let m existing 
facilities be located at known distinct points P][_j*** 5^21 
and. let n new fac.ilities be located at points 
in the plane. Let ?^i^ represent the distance between 
the locations of new facility 3 and existing facility i 
and d(X^,^) be the distance between the locations of the 
new facil.lty j and k. Let the annual cost per uni.t distance 


betV'/een new facility 3 existing facility i be denoted 
by Wj^j with v^^ being the correspond.lng annual cost per 
unit distance between neif facilities 3 and k. The total 
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annual transportation cost associated with new facilities 

located at is given 

n-1 n n m 

= I I V &ix ,x,) I T 

3=1 3=1 1=1 

(4.1) 

because in order to evaluate the annual eost due to item 
movement "bet^Treen new facilities, it is only necessary to 
su'.m over those values of 3 which aro less than Iz and over 
values of k frora 2 to n. 

Somewhat less geometrj.cal insight is available for 
multi- facility location problems as compared to the single 
facility problems* The construction of contour linns is 
no longer possible, because n being at least two, the number 
of variables increases. 

It is the costs of materials handling be tore en 
the new facllitl.es which distinguish the multi- facility 
loccition problem from the single facilit3r problem. This 
follows, hecause when all v.-j^^s are zero, (4.1) majr be 

(4.2) 

(4.3) 


writ'ben as 

n 

f C I f cxp 

■ 3=1' 

m : 

where f, (X..5) — Z d(X. , Pj) 

3 3 3 m 3 V 1 


i=l 


We obseirve that (4,3) just defines a s-iiigle facJ.lity 
total cost expression while (4,2) is the sui-i of n different 
single facility cost expressions. Since the location of 
one new facility has no effect upon the cost of locating 
other new facilities, we have 

n 

min f (X 3 _,..,, X^) = J min f^.(X.) 


That is, least— cost locations of the new facility 
may he found hy solving n single facil3.ty location problems 
independently. Thus, the terras v in (4,1) give the 
multi- fac5.1ity location problem its special character. 

The new facilities j and h are said to have an 
e::change when v^^^ is non-zero, and to have no exchange 
when is zero* It will alwa 3 rs be assumed subsequently 
that each new facility has an exchange with at- least ore 
other new facility. 


I'Tot only shall we assume that there is an exchange 
between new facilities, but also it will be assumed that 
there is an exchange between new and existing facilities. 
As a motivation for this assumption, consider a situation 
5_n which there exists a collection of new facilities that 
have exchanges only among those new facilities within liie 
collection. Where should facilities within the collection 
be located? Obviously, all tie facilities must be located 
at the same place and this place can be anyifherc in the 
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region So there is no location problera. In subsequent 
discussion we shall exclude such meaningless situations. 

Similar to the sn.ngle facility location probloias, 
in this case also, we must tahe and d(X.,P^.) as 

appropriately defined distances. The three types of 
distances discussed previously will be dealt with l.n this 
chapter also in fo iiiiu3_ating the continuous multi- facility?' 
lo cation p robleni s • 

The discrete version of the multifacility ].ocation 
problem has been treated by Cabot, Fra,ncis, and Stary (3), 
hinds che die r and Moore (2), Francis (14) , Mesolowslcy and 
Love (44), and Pritsker and Chare (34) for rectilinear 
distanice, by Eyster and White (S) and VThite (45) for 
Squared Euclidean distance 5 and by Francis and Cabot (15), 
Love (26), and iliehle (30) for Euclidean distance. The 
methods of. Eyster, White, and ITierwille (8) and We so low sky 
and Love (42) can be used for rectilinear problem as we. 11 
as Euclidean problem, 

4.2 C ont inuous Deman d Mu lti-Facility Loca tion Proble m i 

Just as we extended the discrete location problem 
to continuous problem for the case of single Tiew facility, 
we shall consider similar extension of discrete multi- 
facility location problem to continuous multi-facility 
location problem. We sha.ll assume that all the .new^ 
facilities have demands ^arising all over 3n tlie region 

. t ■ hh.. " ' 
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under consf.de ration* For solving the continuous version 
of the multi- facility location problem, ve shall adopt 
•'iTt'J'o approaches, sfmilar to those used for solving continuous 
s.lngle facility problems* Incidentally, the objective 
functions foimula ted by these ti-ro approaches can be minijnized 
using the same methods for all the three types of distances. 
So we shall give the methods of solution for the first 
approach only. 


4,_ 3 Rec tangular Area Aup roac h t 

In this approach we divide the region under 
consideration into several rectangular a,reas having cons- 
tant denand.s. It is to be noted that the region will have 
to be df.vided in as many ways as there are the new factf.f.ties, 
i,e, 5 we shall have to divide the region in n ways, one 
for each facility. 

Lot p. be the total number of re cta.ngular areas 
in which the region ha.s been divided for the demand of jth 
new facility and u., be the constant demand for jth facility 
arising in kth rectangular area, k =1,,,,, p^. Let us 


call this area 

may be written as 
n-1 

j=l 


Then the total transportation cost 


n 


I 


11 


jk ^ 


n ^2 

I I dCXpX) dz dy 

3^1 ija ^3k 


(4,4) 
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ifliere Z is any point in the plane. 


4,3.1 Re ctiline.ar , Disteaice Problem - 

In this 

reduces to 



n~l 

n 


f(x y 3 _}...jyj^) = 1 

3=^1 

1?=3+1 



n 

P3 

1 yj-yji 

+ I 

1 ^ 




A 

Die 


[ |x. - xj + |7^-y|3fe dy (4,5) 

Since the above expression is separable in Xi* s 

tJ 

and y.*s, we can find their values independently. 

D 

(4,5) may be rewritten as 

where 


L^( 5 • • » 


n 

I I W 1/ ix -x| dx dy + 


5^1 lc=l 
n-1 n 


A., 

3k 


3=1 1 ?= 3+1 


(4.6) 
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and 


n p_ 


,y^) 


tj 

^ ^ 1 y.--y! te dy + 


3=1 1^=1 

n-1 n 

i i '' 

3=1 fc=3+l 


A., ■ 3 

3k 


3k • 3 


-II 


(4.7) 


We shall discuss the miniinization of (4.6) orO.jr, 
since the nim:jnizat5-on of (4,7) is an entirely similar 
procedure. 

Let ^ ^^3k’ *^31? boundaries 

of the rectangu3.ar area in x and y directions respec- 
tii;ely. Then (4,6) reduces to - 






n 

'’-3 

I 

I « 

3=1 

fc=l 

n-1 

n 

I 

I ^ 

3=1 

k= 3+1 

n 

I'3 

I 

I u. 

3=1 

k?=l 

n-l 

n 

I 

I 

It 
‘ H 

k=^l 

1:' 

, -i, ■: 




'3k ^^'3k ** *^31^^ ^ ^ ^3"^ ^ *** 

^•1 
3k 


3k ' 3 k 


^3k 


3k 


/ j j ds + 


^3k 


5 k ' 3 


1.1. T. 

fFIstTiitAL IISRA^V 

m i *1 C C 



The minimization of (4*8) can he easily done hy 
the univariate method or some other technique of non- 
linear pro grarami.ng. Since the oh jective function is 
convex, the convergence nil! he towards the global minimum. 

However, if we can tolerate some error, an approxi~ 
nation to (4,8) can he minimized directly h^?" adopting the 
procedure given helow. In this procedure, the terms 

lx. *- I and 1x. - x I are taken as [(x.-x,J^ +e.l and 
2 

[(x. - x) + e- 1 respectively, e being a positive vaJmed 
small constant. 

He put the above mentioned approximations in (4,8), 

and then differentiate it partially with respect to x. ; 

3 

3 = l,,.,,,n*s obtaining 


3f^/8X^ = 


n 


X 

i 

k=l 



1^3 



^3 

f 2 

X 

(x^-x) 

I 

V J 

+ e ] 3 

k=:l 




,n 


When this is equated to zero, for getting a minimum value 
of f^j we have 


n 

l Pk V [ + 

^3 n ^ 




p, 

i 


TJ... / xdx/C(x.-x) + 

3K a 3 

3X 


b^ 


’3k 2 

J dx/ [(x.--x) + e] 
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The above expression can be used as an iterative 
formula for obtaining s. Similar formnla can be derived 
for y s. 

tJ 

4j'iLu^ula^iaxeAJpxjLi.deaii.^D^^^ In this case 

the total cost will be given by 

n-1 n 

,y^^) - I I 

n p . 


l I // [ Cx.-x)^+(y.-y)^] dxdy 


3=1 lc=l 


3=1 te=3tl 
3 ■ "'3 




A., 

3 k 


Since this expression is separable, ve can indepen- 
dently find out the x aid y coordinates of the new faciliti.os, 

X co-ordi.nates can be found out by the minimisation of 

n-1 n 

= I i + 

3=1 15=3+1 

n Va 


I I iiji, // tedy 

3=1 1-1 tk 

For a minimum value of f^, we must have 


(4,10) 


3f^/ 3x^ = 0 ; 3=1,,., ,n 


P. 


n 


r 


I U,,, // 2 (^. - .tady + I 2 V (X 


:0 


3k 


bk ■ 


lc=i 3k J 

^3 


1^1 
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Or, 





n 


A,, 

1?=1 Dk 


// X fe dy + I V.. X 
^ 3 k k 


"3 


lc=l 


n 


11 


3 k 


// dx dy + V. 


A, , 
k 3 


k=l 


3 'k 


i ~ 1 , ♦ • * jii» 


Similar expressions for y.^s can "be deteruiined by 
using tlie same arguments. Hence we observe that we shall 
have to solve ti-^o systems of simultaneous linear equations, 
each system having n equations. 


443,0 Suclidean Dista nce Probl.em - The total transporta- 
tion cost of the Euclidean distance problem is 

n~l n 2 

f(xq,...,x^| yq,..*,yj^) = I I V^^ [(x^-XjP + 

3=1 is?= 3 +k 

/ -2 V 2 

^ Pj 

' ^ ■' 2 ' 'i 

+ I I ki V // [(sr.-x) +(y^-y)^ ] <lxdy 

3 k A 3 J 

3=1 l 5 =l jk 


We observe that our objective function is a non- 
linear convex function of the independent variables, x. 

3 

and y . , 3 = 1 , • . • ,n, so noivlire ar p ro gramming technique 3 
can be successfully used to obtain a convergent global 
optimal solution to the problem. 


"i. - 1 I 


‘■0.- 
'C 




..aP: 
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Al'csriia.tely w© ca.ii solv© "bliis t)i“o'b3-eiji in ili© 
j.o3.1oi'iiiig "Way also* He Imow that for a iiiinjjinitii value of 
f j we must have 


3f ^ 3f 
3 71 




= 0, j = 1 


>• * » } 


n 


P. 


or 


I tk II 


li=l 






{ (x.~x)^+(y2-y)^] 


-TT- dy + 


n 


I 


3 k 


( Xj~Xk^ 


or X. 


V 'k -^1-^ 

^1. Ilk .. 1-i . llr 


^3 


n 


I // + 1 

3k A, “ 


3k 


k=l ^ li:=l ^3k 

n 


Siiiij.larly y 


3 


--tel 




^ . ^iktefc 

-7-<3:^fe - + I- -3~ . . 


tel B., 
tel 3k 


n 


I u., // •4saz+ I 

tel -^jk 3 


V 


•where 3. = [(x^-x)^ + (y-.-y)^ 1 

■■ 2 ' ' ' 2 


•I 


tel 

te3 


^3k 


and 


1 

2 2 




(4*11) 


(4*12) 
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If, for any 3 and k, (x^ ,y^) = (x^^, y^ , then 
iDecoiiies zero and ( 4, 11 ) and (4,12) are undefined. 
Consequently, ve see that difficulties arise when the 
locations for any two new facilities coincide (matheiiiati- 
cal3.y) with each other. If there were some guarantee that 
the optniumi locations of any ti-ro new facilities would 
never be the sajae, then (4.11) and (4,12) would still give 
the solution to our problem, because these equations, then, 
could be solved iteratively to give the optjmal loce-tions 
of new facilities. Unfortunately, there is no such guarantee 
a-vailable. Consequently, a modification of this approach 
is required. The proposed modification is to use an appro~ 
x5mation of instead of D-j, in (4.11) and (4,12). This 
approximation is to be such that even if for some 3 and ]c, 
(x., y^) = (x^^, YiP? does not equ.al zero. The proposed 
approximation is to use 

1 

p 



in place of where e is an arbitrarily small, positii?B- 

valued constant. The valiae of e can be de'teimined by 
a compromise beti-reen how quiclcly we want to get the 
solution of our problem and the accuracy of that solution. 


The final iteration equations can be given by 







p * 

D 


n 


:x 


(M) 


I "3k // I 

__ " , Iq^j 

I u // mz + I 
l.=l -^3kS/« k=l 


■“Tnrr 


(h) 




^jk 


3k 


P 


n 


D 


^3 


(h+i)^ S: 


I ff- I .J}1^ 

_^__f_3_ 3k 

Pi n „ 

^ * I 

k=i 3 k 


(h) 


I ’^jk //,.(k) 

^3k 3 


151=1 


1q^3 


■■ De n si ty Func tion A]3p,rpaj;1:i j 
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In this approach we define the demand density 
functions f .(x,y), 3=^1, ,.,n, over the entire region under 

tJ 

CO ns.1 deration for all the new facilities to he located* 

Then the total cost of transportation will he given hy 
n-1 n ^ 


f(Xj_,..*,2^) = I 

3=1 


1 //fj(x,y) 

1^3+1 ^ 


d(X^,X)(ix:dy 

where A is the area under consideration! and d(X^,Xj^5 

and d(X.,X) denote appropriately defined distances, 

, .< 3 ' ■ 

Bie total co st 

function,; in this casej becomes 


n 
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n~l 


yi,...,yn) = I r V [ I X -x } + 
3=1 1^3+1 





n 


1 yj-y!: 1 

] + 

I 




3=1 

uhi.ch 

is sepa: 

rable 

In 

CO rre 

spondi.ng 

to X 




n- 

i 

f^( x^ 

- ) 

J ♦ • * 5-''22'^ 

ft 



i ' -3'"’'^' I nj I ■ ' ^3- 


Idxdy 


n n 

^ ^3k I I ^ 

3=1 1~3+1 ^ 


f .(x,y) lx— x| dx dy 


wlij.ch can lie either m3.rd.mized by the univariate laethod of 
nonlinear progranmiing or by the procedure given in section 

4*3,1. 

4,4.2 Squared Bluclidean Dista n ce .PrpbJ_.erii « In the case 

of squared Euclidean distance the ejnression for total 

transportation cost vill be given by 

n-1 n 

r 2 

fC^i j » • « y3 _5,**jyj^) - I 1 ' [(x— Xj,) + 

3=1 ic=3+i 
o ' n 

2 2 

(y^-yJ ] + I // f^Cx,y) l(x— x) + 

3-1 ^ 

(y^-y)^]<3xiy 

Tills function can be minimized in a manner similar to one 
given in section 4,3,2, 
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“ The objective function 
for the Buc3.idean distance problem is given by 

n-1 n 


fCx- 






j=l l^j+1 


hi/ 


1 

, 2 . ^ 


n 


(7^-7^) ] + I // f^(x,7)[ (x^-2:) +Cy^-3r) 3 

3=1 ^ 


d2cdv 


This function is also convex and amenable to solve by 
a.ny non-linear programming technique which could bo used 
for unconstra.lne d problems* But a c?:-ose a.pprood.mat.ion to 
the above objective function can be minjmized by a procedure 
similar to that of 4.3,3, 



Chapter 5 

The problem' which we shall deal in this chapter 
is-aspecial type of discrete demand single facility location 
problem,,, rather thana continuous demand location problem. 
h‘e, still, propose to discuss it here because 
it is very much sjmilar to the cont^mious demand 
single facility location problem as far as its nature and 
method of solution are conce3?.!ed. Piirtherrnore , it is an 
important problem wh.i.ch arises in real life quite frequently 
a.nd is an extension of the discrete demand single facility 
location problem. 

All the previous research work on the discrete 
vdomand facility location problems asstmies that the locations 
of the destinations or the existing f ac.ilitie s ' are deter- 
ministic. In the present formulation, the destinations are 
no longer predeteriiiined pointsj but random va.riables with 
given probability distributions. The problem we now wish 
i:.o solve is to find the location of a point, the sum of 
the expected values of whose weighted distances from all 
destination points is a minimuffl. This model i.s useful in 
all real life situations where 'some uncertainty is 
associated with t|ie locanihii of the destinations. 



In order to formulate the prohlem, let ns assume 
that I'/e are given n fixed destinations 5 whose 

co-ordinates are (x., yj, 3 = l,...,n, and each P . has 

associated with it a probability density function f.-Cx.- ,y^.) . 

J J 

The demand at jth dest?_nation, i.e., at P,. is assumed to be 
w . and the new facility is located at a point X irhose 

t) 

co-ordinates a,re (xj y)» Then the e3:p)ected va.lue of the 
weighted distance between new facility and 3th destineution, 
i.e,, P. will be given bjr 

J 

'iW, Pj) -i'jCxj, yp tej dyj 

and hence our objectj.ve function considering all destina- 
tions will be 
n 

00 ^ 

f(x) t= J !^J_^ yp ^3* 

3=1 

which is to be minimized with respect to x and y, 

We observe that the foim of f(x, y) is similar to 
the ob jective function for conttjnuous demand single ’facfiity 
location problem formulated by the density function approach 
with the only exception that in ( 5 , 1 ) one additional 
siimiiiation sign appears in the beginning. But, fortunate 3 .y, 
this difference does not alter the situation much and the 
solution procedures used for continuous demand single 
facility problem can'be adopted with slight modifications. 



Th6 modifications for all tlie three distances are presented 
in the following paragraphs* 

1* ^ When rectilinear 
distances are used (5.1) tales the form 

f(Z) = fi (x) + fg (y) 

. . n 

- ■ 

where f-j^Cx) = I ] |x>-x| T.(x,-) dx. 

3=1 

n 

f2(y) = 1 wj LVj-yl 
3=1 

and FjCxp and G^Cy^') are marginal density functions of 
f^Cxj, y^) with respect to x^ and y^ respectively* The 
optimal values of x and y can he ohtaiiied hy the indepen- 
dent minimization of f 2 .(x) and fgCy). 

By arguments similar to those used to solve 
continuous demand si.ngle facility location problem by 
density function approach, we obtain fol.lowiiig condition 
of optima.lity for x* Similar condition can be had for 
y also, 

11 . 

Y wYf =0 ^ ^5,^ 

/„ 0 I- «J ■ ^ 

i=i 

where H^(x^) ^jF^(x^)dx^, 3 ® l,,»,jn 

,2)!i-S an algebric equation' in" x and can be solved (at least 
numerically) to yielsl;^^ jOPt value of It. 



2* S^axeA JlucjLldean.PxoMem. Proceeding exactly on the 
same lines as we did .In Chapter 2 with density function 
approach, we can ohtaDji the fo3J.owing optimal values of x 
and y. 


X = 


and y 


n 

n 

■<r~ 

„ • 

f— W4 m / 

2 x'A ^ 

3=1 ^ 

Z- 


n 

n 

^ , ''3 “y/ 

3=1 


H 

11 


w. 

J 


w. 

;) 


where rn^ _ and m are respectively the means of the margi- 
nal density functions (xj) and G^Cy^, This shows that 
the prohahilistic problem for squared Euclidean distance 
can be solved by assuming it to be a deteiministic one, 
with the co-ordinates of the existing fac5J.j.ties tahen as 
the mean values of the probability density functions 
associated with each destination* 


iE!upl^.d.e,an ,Djntanqe„, Prpbl_em« When Euclidean distances 
apply, we can devise an iterative scheme similar to the 
one used for continuous d.emand single f anility location 
problem,. The iteration formulae will be 


(lw-1) 


n 

z: 

3=1 


OO 


■ oa 


3 „ "'3 ayj/ dj 


(h) 
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silngle and nmltifacri.lity location problems for all the 
three tjrpes of di. stances have been treated by this approach. 
It is to be pointed out that this a.pproach is easier and more 
acciirate because dividing a region into rectangular areas 
uith some demand is a difficult task, Involves quite a 
sufficient amount of error, and no systematic procedure is 
available to do it. On the other hand, the detemiination 
of a demand density function can be done by ev0.1uating the 
demand at points of a grid constmcted over the region, and 
then by interpolating beti^een them, which is a systema,tic 
procedure easily done by the techniques of numerical 
analysis. The accuracy of this procedure can. also be 
fudged. Furthermore, the solution procedures bf highly 
similar na-ture are obtained in both of the two approaches, 
especially for squared Su cl idean and Buclidean problems, 
hence it is better to use density funct3.on approach, 

1 nrobab ill Stic model . of , tho discrete demand 


single facility location problem has* also been developed. 


In this 


model the locations of destinations or existing 


facilities are random variables wii^, known probability 
distribution functions. It turns out that this problem can 
be solved exactly in the same.^s^p^r.as Idie continuous 


demand single .facility, 'solved by the 
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The procedure^ for developing the contour 15_nes 
for these prohahilistic models has not been discussed but 
it is hoped that this can be done easily* Itirther the 
methodology adopted for the probabilistic model- of the 
discrete location problem can be extended to the multi- 
fa,cility location problems. 

Some continuous demand facility location problems 
arising in practice may have some . constraints upon the 
location of the new facilities. If these constraints 
exclude some part of the region, str'd.1 mailing some conti- 
nuous space available for locating the new facilities, then 
the problems can be formulated as in this thesis (of course, 
they will have to be solved by non-linear programming; the 
special methods devised will not work), but if the constraints 
are such that they leave only some discrete points available 
for locat5.ng the new facilities, then new methods will be 
required to handle the problem. Therefore, further 
research efforts need to be directed towards the exten?- 
sion of the discrete spa-ce location problems discussed jn 
fS-rst chapter for continuous demand situations. 
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Another scope for futiire research is to extend 
tjie continuons demand multifacilitj;" location problem to 
one involving sUnilar facilities. Both nncapacitated and 
capacitated versions of the shnilar facility problem may 
be undertaken. One can also extend the models presented 
in this work to those which involve both discrete and 
continuous (point and area) destinations, or to those 
Involving a combination of various distances because some 
facilities may be connected by a straight line distance 
whereas others, by rectangular aisles. The location theory 
models developed, in this thesis pertain to two dimensional 
space. Attempts shou.ld be made to extend these models to 
three dimensions to represent a situation in which facili- 
t.les are to be loca.ted in a multi- stoi’eyed building at 
different floors. 
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